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In this paper we calculate C(m*) corrections from holonomies in the Loop Quantum Gravity, 
usually not taken into account. Allowance of the corrections of this kind is equivalent with the choice 
... of the new quatization scheme. Quantization ambiguities in the Loop Quantum Cosmology allow 

00 ' for this additional freedom and presented corrections are consistent with the standard approach. We 

\ apply these corrections to the flat FRW cosmological model and calculate the modified Friedmann 

. equation. We show that the bounce appears in the models with the standard 0{^'^) quantization 

(N . scheme is shifted to the higher energies pbouncc = 3pc- Also a pole in the Hubble parameter appears 

\ for Ppoie = f Pc corresponding to hyper-inflation/ deflation phases. This pole represents a curvature 

, singularity at which the scale factor is finite. In this scenario the singularity and bounce co-exist. 

' Moreover we find that an ordinary bouncing solution appears only when quantum corrections in the 

lowest order are considered. Higher order corrections can lead to the nonperturbative effects. 



O^, I. INTRODUCTION 

. 

Strength of the gauge field F in some point x can be obtained form holonomy calculated around this point and taking 
limit of the zero length of the loop. Loop Quantum Gravity (LQG) is kind of gauge theory describing gravitational 
ps) . degrees of freedom in terms of gauge field A which is elements of su(2) algebra and conjugated variable E which is 
^ ' elements of su(2)* algebra To quantise this theory in a background independent way one introduces holonomies 
CO , of the Ashtekar connection A 

o 



hc[A] = "Pexp / A where 1-form A = r^A'^dx"^ 

J a 



(1) 



where = — fcr; (ci are Pauli matrices) and conjugated fluxes 

oo ■ 

o 



Fi;[E]^ dF' where 2-form dF, ^ eabcEyx'' A dx" (2) 
Js 

■ I 

■ as new fundamental variables 0, Q . Other variables like the field strength F should be expressed in term of these 
elementary variables. As we mentioned at the beginning the field strength can be expressed in term of holonomies. 
However, another aspect of loop quantisation starts to be important here. Namely, an area operator possesses a 
discrete spectrum with minimal nonzero eigenvalue A So we cannot simply shrink to zero the area enclosed by 
loop. Instead of this we must stop shrinking loop for a minimal value corresponding to the area gap A. This effect 
leads to quantum gravitational corrections to the expression for classical field strength. The expression for the field 
strength as a function of holonomies have a form 



tr 

= lim <' -2 



(3) 



where the limit ^ ^ Jl corresponds to the minimal value of the area gap A. However, this formula is adequate only 
when ©(/i^) terms can be neglected, i.e., in the classical limit. In fact these terms, which form infinite series, are a 
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function of F and A. The expression for the F as a function of holonomies should be therefore obtained by solving 
this equation in terms of the first factor on the right side. In the classical limit /i ^ terms 0{fj,'^)/fi^ vanish and we 
recover a classical expression for the field strength. Until now in literature the first order quantum correction to field 
strength has been investigated. It means that terms C(/i'*) have been neglected. This approach was dictated be the 
choice of the simplest quantization scheme. Namely, as it has been shown by Bojowald [a, Hill, the precise effective 
Hamiltonian must be a periodic function of the canonical variable c. The simplest form of this function we obtain 
when we perform the regularisation of the expression for the classical field strength cutting off the terms 0(/i*). This 
is a standard procedure in the Loop Quantum Cosmology. 

In this paper we calculate and study another non-vanishing contribution what is in fact a choice of the different 
regularisation of the expression for the field strength. It means that we hold C(^^) factor, which is a function of F, 
and we solve equations for as a function of the holonomies. This approach is equivalent to the choice of the new 
quantization scheme what is allowed due to quantization ambiguities. 

The organisation of the text is the following. In section |TT] we calculate expression for F as a function of holonomies 
in 0{fj,'^) order. Then in section ITTTl we apply this result to the flat Friedmann-Robertson- Walker (FRW) cosmological 
model. We show that obtained correction have important influence for this model. In section |V] we summarise the 
results. Finally in the Appendix we give some basics of Loop Quantum Cosmology connected with the subject of this 
paper and explain the employed notation. 



II. HOLONOMY CORRECTIONS 



From the definition ([T]) we can calculate holonomy for homogeneous model in the particular direction °efda and 



the length ^iVq 



1/3 



Icos(^) +2nsin(^) 



(4) 



From such particular holonomies we can construct a holonomy along the closed curve a = Dij as schematically 
presented in the diagram below 




o an 



and can be written as 



(5) 



where we have introduced 



(6) 



Factors Bi are elements of 5u(2) algebra so to perform product of exponents in equation ([5]) we need to use the 
Baker-Campbell-Hausdorff formula 



e-^e^- = exp ^ X + r + y] + 1 ([X, [X, Y\\ + [y, [r, x\\) - [x, \x, y\\\ 



(7) 
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To calculate ©(/i^) correction the elements of the expansion written above are sufficient. Applying this formula to 
equation ([5]) we obtain 

r 3 4 4 ^ 

hl^l - exp U'[B,,B,] + yfS. + B„ [B^B,]] - ^[B„ [B,, [B,, B,]]] + ^[B^ + B,, [B, + B„ [S„ S,]]] + 0{ti')\ 



= I + fi'[B„ B^] + ^[B, + Bj, [B„ Bj]] ~ i^[Bj, [B„ [B„ B^]]] + ^[B, + B^, [B, + B^, [B„B. 



12 
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-[B,,Bj][B,,B,] + 0{fi^). 



(8) 



Now, multiplying this expression by Tk , using definition ([6|) and taking a trace of both sides we obtain 



tr 



(9) 



We mention that {ijk} are external indices and the Einstein summation convention is not fulfilled. The introduced 
parameter Ch corresponds to the effective canonical variable c which is expressed as a function of holonomies. With 
use of equation ([4]) we can directly calculate the left side of equation ^ , we obtain 



tr 



-^sin^ (/ic). 



Then, using properties of Ti matrices we obtain 

i^^'ch - /i^Ch + sin^ (/ic) = 0. 
The 0{^^) order contribution simply vanishes. The solutions of this equation have a form 



1± 



^l-|sin2 (^c) 



-U2 



When we expand the square in the solution for cf-^_ we obtain 



sin(/ic) 



1 sm^ (fic) 
3 /i2 



(10) 



(11) 



(12) 



(13) 



The first factor of the expansion corresponds to the known case when 0(p'^) corrections are ignored. We can easily 
check than the classical limit /i — > 0, ^ is recovered only in the case. The case c^^ should be therefore 
treated as unphysical. However, as we will see in the next section, both solutions lead to the same modified Friedmann 
equation. So we can keep both solutions. 

Finally the expression for the effective field strength has a form 



1± 



We have performed here the limit /i — > /i where 



Jl-|sin2 (/2c) . 

V 2 o to J 



2-2y2/3 



(14) 



(15) 



For details of this limit see papers 0] or appendices to the papers d, [T^, [Tl| . 

As we have mentioned earlier the precise effective Hamiltonian must be a periodic function of c. In our case the 
effective Hamiltonian has a form H^s \/lp[c^_ where the Cy^_ can be expressed as 



(2n)! 



2u2 (2n - l)n!23"-i(2i)2" 



[exp(i/^c) — exp(— i/ic)] 



2n 



(16) 
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As we see this function is periodic, and forms an infinite series numerated by integers. However, this infinity is allowed 
in the frames of Loop Quantum Cosmology. The obtained effective Hamiltonian is correct however is not given by the 
simple function as we should expect for fundamental expressions. However, we should to keep in mind that we are 
looking for the effective Hamiltonian and there is no circumstances that such a Hamiltonian must have mathematically 
simple allowed form. 

In the next section we will use the calculated effective field strength F for the FRW k = cosmological model. 



III. APPLICATION TO FRW fc = 



With use of equation (fTl|) we can derive the effective Hamiltonian for the flat FRW model in the form 



1± 



Hcs — 



(17) 



For details we send to the appendix. This Hamiltonian fulfils the so called Hamiltonian constraint H^s — 0. From 
Hamilton equations we can calculate evolution of the canonical variable p 



and with use of (fT7|) we obtain 



P=T 



dc 



\p\ 2sin(^c) cos(^c) 



1 — I sin^ (/ic) 



(18) 



(19) 



Applying equation (fT9|) . the Hamiltonian constraint i/cff = and definition of the Hubble parameter H = we 
finally derive the modified Friedmann equation 



ff2 



_P_ 

3pc 



(' - if) 



where we have introduced 



V3 



(20) 



(21) 



PI 



As we see obtained equation does not depend on the sign ± in the Hamiltonian. An analogous equation in the lowest 
order has been calculated earlier fsj and have a form 



fr2 



SttG 



1 



P_ 

Pc 



(22) 



This equation lead to the bounce for p — p^. An analogous bounce is also present in the derived model (|20p. however 
now the bounce is shifted to the higher energy densities 



Phc 



3/Oc 



Another important property is the pole in the Hubble parameter for 

3 



Ppolc 



(23) 



(24) 



as we see from equation (1^01) . We show these features in Fig. [T] In the left panel we present as a function of 
energy density for 0{p^) and 0{p^) cases. In the right panel we compare evolution of the Hubble parameter as a 
function of p for the radiation dominated Universe (p cx Ijp^). 

In both case we observe a nonperturbative feature, namely the pole for p — |pc- This fact indicates that higher 
order corrections from holonomies can have important influence for dynamical behaviour for small values of p. It is 
clear from a parameter of expansion which grows for small values of p. For large p the classical case is clearly 
recovered, however behaviour for small values of p is highly complicated. Namely, as our study suggests higher 
order terms of expansion have nonperturbative influence for dynamics and this fact can seriously complicate a simple 
bouncing universe picture. We investigate this issue in the next section. 
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FIG. 1: Left : Evolution of H^^^2^ (red, bottom curve) and ^^0(^,4) (green, top curve) as a function of energy density p. The 
parts below zero value of are unphysical. The p is in the units of pc. Right : Evolution of (red, bottom curve) 

and H^i^^i^ (green, top curve) as a function of p for the model with radiation (p oc 1/p^). On the both panels a physically 
admissible region corresponds to > 0. 



IV. QUALITATIVE ANALYSIS OF DYNAMICS 

The advantage of qualitative methods of analysis of differential equations [l^l is that we obtain all evolutional paths 
for all admissible initial conditions. In this approach the evolution of the system is represented by trajectories in the 
phase space and asymptotic states by critical points. We demonstrate that dynamics of the model can be reduced to 
2-dimensional autonomous dynamical system. These methods allow to distinguish a generic evolutional scenario. 

In Fig. [2] we show the phase portrait for all admissible initial conditions (all values of total energy E of the fictitious 
particle moving in a 1-dimensional potential proportional to p^) for the model with the free scalar field. 



P' 




FIG. 2: The phase portrait for all admissible initial conditions. The blue line is representing points at which trajectories pass 
horizontally through the inflection point {hyper-mflationary /deflationary phases). 

The physical trajectories are situated in the region at which E — V is non-negative or p is larger than some 
minimal value and p' is zero. The physical trajectories lie in the non-shaded region bounded by a zero velocity curve 
which represents a homoclinic orbit. Of course the whole system is symmetric with respect to the reflection {H is 
changed in ~H). The blue line on the phase portrait represents points at which trajectories pass horizontally through 
the inflection point {hyper-inflationary /deflationary phases). Therefore, evolution comprises the bounce solution 
interpolating static phases of evolution, see Fig. [3l There is the intermediate phase of evolution at which we have 
the inflection point at the diagram of p{t). Note that it is the singularity state with rapid growth of the scale factor, 
we call this phase the hyper-inflation. It is important to note that energy density is finite during transition through 
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these singularities. Similar finite scale factor singularities has been studied recently by Cannata et al. [l 




rgy density 



scale factor 



^ time 

classical 
singularity 

FIG. 3: The schematic picture of the evolution of the model. The top curve (green) represents the energy density of the 
free scalar field. It is worth to note that this energy density is finite during whole evolution, even during transitions through 
singularities. The bottom curve (blue) presents schematic evolution of the scale factor for the investigated model. The dashed 
curve represents the classical evolution which is not realised in the presented quantum model. 



In the phase diagram in Fig. Owe adjoint a circle at infinity in standard way via Poincare construction. Note that 
all trajectories are starting from unstable node - representing static Einstein Universe and landing at stable node 
representing another static Einstein Universe. 

Parisi et al. pointed out that stability of Einstein static models in high-energy modifications of General 

Relativity is important from the point of view of so-called Emergent Universe scenario [l^. Mulryne et al. [T6| 
investigated the stability of the Einstein static model and they found that LQC Einstein static model is representing 
a centre type of critical point on the phase portrait. As it is well known such a critical point is structurally unstable. 
Note that in our case this static universe is representing a node type of critical point. This modification of stability in 
presented model has important consequences for the Emergent Universe scenario, since as it is well known in General 
Relativity, a static Universe is unstable and is represented by a saddle type of critical point and therefore it requires 
the fine-tuning. Moreover corrections considered lead to a regularization of the big-bang singularity [13 ■ 



V. SUMMARY 



In this paper we have calculated another non-vanishing contribution to the quantum holonomy correction in Loop 
Quantum Gravity. Quantum correction of this kind appears when we express the Ashtekar connection A and field 
strength F in terms of holonomies. The source of quantum modification to classical expressions is a non- vanishing 
area enclosed by loop as the result of existence of the area gap A. 

We had applied obtained corrections to the flat FRW cosmological model and the we have calculated resulting 
quantum gravitational modifications to the Friedmann equation. The holonomy correction in the lowest order to the 
flat FRW model has been calculated earlier 0, [l^, [l^ and extensively studied [13, 13 ■ These investigations uncovered 
existence of the bounce for energy scales pc- In this picture the standard Big Bang singularity is replaced by the 
non-singular Big Bounce. Calculations performed in the present paper indicate that the holonomy correction in the 
next non-vanishing order holds this picture. Namely, the initial singularity is still preserved. However the bounce 
appears now for higher energy density pbouncc = Spc- Another important new feature is the appearance of a pole in 
the Hubble parameter for ppoic = fpc corresponding to hyper-inflationary/deflationary phases. This leads to more 
complicated dynamical behaviour at these energy scales. 

We showed that the generic evolutional scenario for the model with the free scalar fleld starts from the static 
Einstein universe then recolapse passing through the curvature singularity with a flnite scale factor (hyper-deflation) 
towards the bounce and goes in the expanding phase through the second curvature singularity (hyper-infiation) and 
ends in the static Einstein universe. During the transition through the singularities the universal critical behaviour 
H (X \p — ppoiel"^ holds. Therefore in the presented scenario the bounce connects these two finite scale factor 
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singularities. 

As we see, the higher order quantum correction in LQG can have an important influence on dynamical behaviour of 
cosmological models. It is not unlikely that the non-singular bounce appeared in the lowest order can be only an artifact 
of simplifications and can disappear when whole contribution will be taken into account. Further investigations of 
quantum corrections from LQG arc still necessary. We conclude that higher order holonomy corrections and resulting 
different quantization schemes should be also seriously taken into account in considerations. 
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APPENDIX A: FLAT FRW MODEL IN LOOP QUANTUM GRAVITY 



The FRW fc = spacetime metric can be written as 

ds^ = ~N^{x)dt^ + qabdx°-dx'' (Al) 
where N{x) is the lapse function and the spatial part of the metric is expressed as 

qat = = a'itrq,, = a2(<)%°c.r^^ (A2) 

In this expression is fiducial metric and "w* are co-triads dual to the triads °ef, °Ci;'(°ej) = (5* where °cj' — "uj^^dx"" 
and °ei = °efda- From these triads we construct the Ashtekar variables 

a: r^ + 7/C-cv'/'X, (A3) 

Ef = v^d^e^=pV'^'V^X (A4) 

where 

\P\ = (A5) 

c = -faVfl^^. (A6) 

Note that the Gaussian constraint implies that p —p leads to the same physical results. The factor 7 is called the 
Barbero-Immirzi parameter, 7 — In 2/(7r-\/3). In the definition (|A3|1 the spin connection is defined as 

Tl = -e^^"^-e5(a[,e^] + ie^.ei9[e4]) (A7) 



and the extrinsic curvature is defined as 



Kab = ^ [qab ~ 2D^,Nf,)] (AS) 



what corresponds to AT* := Kabe\. 

The scalar constraint, in the Ashtekar variables, has the form 



1 /• _ EfE'' 



where field strength is expressed as 

F^, = daAl - dbAl + e%A\j^,. (AlO) 
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With use of (fO|l . (|X4)) and (jAlOp the Hamiltonian (M\ assumes the form 




(All) 



where we have assumed a gauge of iV(a;) = 1. Quantum corrections to this Hamiltonian come when we express ^/\p\ 
and in terms of background independent variables. In this paper we had concentrated on the corrections to the 
factor c^j called holonomy corrections. For a short review of quantum corrections we send to the appendix in the 
paper pTlj . 
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